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Summary 


An  analysis  method  for  creating  three-dimensional  meteorological  data,  developed  in  a  previous 
study  (Henmi,  2003),  has  been  adapted  for  use  by  the  U.S.  Army  Integrated  Meteorological 
System  for  the  purpose  of  nowcasting.  This  method  uses  data  regularly  obtained  from  the  U.S. 
Air  Force  Weather  Agency,  which  includes  numerical  forecast  data  generated  by  the 
Pennsylvania  State  University/National  Center  for  Atmospheric  Research  Mesoscale  Model 
Version  5  (MM5)  (Grell  et  al.,  1994),  along  with  routine  surface  and  upper-air  (radiosonde) 
observations. 

This  report  describes  the  method  in  detail,  and  then  shows  how  it  was  applied  to  a  region  in  the 
southeastern  United  States  in  order  to  demonstrate  its  general  utility.  Furthermore,  the  method 
(along  with  two  others)  was  applied  to  a  region  of  Utah  and  evaluated  statistically  using  the  local 
mesonet  meteorological  data.  The  three  methods  used  to  generate  the  statistical  analyses  over 
Utah  were: 

a.  A  nowcast  analysis  using  the  present  method,  and  using  15-km  resolution  MM5  output 
(nest  2  of  a  triple-nested  configuration)  as  background  data 

b.  A  pure  forecast  produced  by  the  MM5,  nest  3  (5-km  grid  resolution) 

c.  A  nowcast  analysis  using  the  Battlescale  Forecast  Model  nudging  method  (Henmi  and 
Dumais,  1998)  and  15-km  resolution  MM5  forecast  output  as  background  data. 

For  the  Utah  domain,  each  method  was  statistically  compared  for  both  a  summer  (2003)  and  a 
winter  (2004)  period.  For  each  method,  a  final  analysis  of  5-km  resolution  was  generated. 

Both  of  the  nowcasting  methods  (a  and  c)  produced  much  better  statistical  results  for 
temperature,  dew-point  temperature,  and  wind  vectors  for  the  5-km  resolution  Utah  domain,  as 
compared  to  pure  forecast  method  (b),  produced  by  the  MM5  5-km  resolution  nest.  However,  all 
three  methods  showed  limited  capabilities  in  simulating  the  high-resolution  wind  fields. 
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1.  Introduction 


The  U.S.  Army’s  Integrated  Meteorological  System  (IMETS)  regularly  receives  gridded  binary 
data  (GRIB)  of  forecast  results  from  the  U.S.  Air  Force  Weather  Agency  (AFWA),  computed  by 
an  AFWA  operational  version  of  the  Pennsylvania  State  University/National  Center  for 
Atmospheric  Research  (NCAR)  Mesoscale  Model  Version  5  (MM5).  The  forecast  data  are 
produced  for  model  domains  of  both  15-  and  45-km  grid  resolution  (and  occasionally  5-km 
resolution),  and  are  generally  reliable  in  terms  of  forecast  accuracy.  However,  the  forecasts  can 
sometimes  evolve  far  from  reality  (Henmi,  2002  and  2003). 

For  Army  operations,  accurate  meteorological  information,  obtained  by  making  maximum  use  of 
currently  observed  local  data,  is  more  desirable  than  information  that  relies  on  fields  obtained  by 
a  previously  executed  numerical  forecast  model.  Nowcasting  is  a  method  of  diagnosing  the 
current  weather  situation.  In  this  study,  nowcasting  is  further  defined  as  an  objective  analysis  of 
meteorological  parameters  that  have  been  currently  observed.  If  the  current  weather  situation 
can  be  accurately  diagnosed,  it  may  be  possible  to  accurately  forecast  or  extrapolate  near-future 
weather  situations. 

In  a  previous  study  (Henmi,  2003),  a  nowcasting  method  that  allows  for  the  rapid  updating  of 
three-dimensional  meteorological  fields  was  developed  and  discussed.  In  this  method,  surface 
and  upper-air  meteorological  fields  are  analyzed,  both  independently  and  univariately.  In  other 
words,  the  surface  observations  and  the  corresponding  MM5  surface  forecast  data  are  composed 
using  the  successive  correction  method  (Sashegyi  and  Madala,  1994),  and  the  upper-air 
radiosonde  data  and  the  MM5  sigma-p  level  forecast  data  are  similarly  composed  using  the  same 
scheme.  After  these  two  steps,  the  similarity  theory  formulae,  MARIAH  (Rachele  et  ah,  1995), 
are  used  to  combine  the  surface  data  analysis  fields  with  the  upper-air  data  analysis  fields.  The 
analyses  must  be  interpolated  both  horizontally  (onto  the  analysis  i,  j  space)  and  vertically  (onto 
the  analysis  terrain-following  sigma-z  levels). 

The  purpose  of  the  present  study  is  to  apply  and  adapt  the  aforementioned  nowcasting  method  to 
a  pseudo-IMETS  operational  environment  and  to  statistically  evaluate  the  method  by  using  Utah 
MesoWest  data  (Horel  et  al,  2002).  Section  2  describes  the  nowcasting  method  and  section  3 
presents  an  example  of  the  method  application.  Section  4  outlines  the  results  of  a  statistical 
evaluation  of  the  nowcasting  method  for  a  domain  over  Utah,  comparing  it  to  results  from  two 
other  methods.  The  concluding  remarks  follow  in  section  5. 
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2.  Nowcast  Analysis 


2.1  Data  Used 

MM5  forecast  data,  along  with  surface  and  upper-air  meteorological  observations,  are  received 
and  added  into  a  database  in  the  current,  operational  IMETS.  The  nowcast  analysis  method 
discussed  here  is  designed  to  perform  computation  using  these  types  of  data.  The  data  that  were 
used  in  this  study  are  discussed  below. 

Data  obtained  for  1200  universal  time  coordinated  (UTC)  on  April  12,  2004  (in  the  southeastern 
United  States)  are  used  in  the  following  examples. 


2.1.1  MM5  Data 

Using  GRIB  MM5  forecast  fields  sent  from  AFWA,  vertical  distributions  of  forecast  data  (per 
grid  point  and  native  sigma-P  level)  are  reduced,  as  shown  in  figure  1 . 
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Figure  1.  A  portion  of  MM5  data  reduced  from  GRIB  data. 

NOTE:  The  top  two  lines  are  the  latitude  and  longitude,  respectively, 
of  grid  point.  Columns  of  data  (from  left  to  right):  pressure 
(mbar),  height  above  ground  (m),  height  above  sea  level  (m), 
temperature  (°K),  dew-point  depression  (°K),  x-component  of 
wind  vector  (m/sec),  and  y-component  of  wind  vector  (m/sec) 
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For  this  study,  12-h  MM5  forecasts  were  used  to  produce  nowcast  analyses  and  make 
comparisons. 

2.1.2  Surface  Data 

Surface  meteorological  data  is  contained  in  the  ASCII  text  fde  named  sf  cmet .  d.  The  file 
contains  observational  data  of  pressure,  temperature,  dew-point  temperature,  wind  speed,  and 
wind  direction.  Figure  2  shows  a  portion  of  sf  cmet .  d  showing  only  temperature  data. 
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Figure  2.  A  portion  of  the  surface  observation  file,  sf  cmet .  d, 
showing  temperature  only. 

NOTE:  Columns  of  data  (from  left  to  right):  number  of  temperature 
observations,  latitude,  longitude,  date,  hour,  and  temperature 
(°C)  per  observation. 
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2.1.3  Upper-Air  Radiosonde  Data 

Upper-air  data  is  contained  in  the  file  named  new  sound,  as  shown  in  figure  3.  Information 
about  latitude,  longitude,  height  above  ground,  date,  hour  of  observation,  and  number  of  levels 
are  given  in  the  first  three  lines.  The  columns  contain  data  on  pressure  (mbar),  height  above  sea 
level  (m),  temperature  (°C),  dew-point  temperature  (°C),  wind  speed  (knots),  and  wind  direction 
(deg).  Missing  data  is  expressed  as  “-999. 0.”  The  file  contains  all  radiosonde  sounding  data 
obtained  from  an  area  of  1,600  by  1,600  km  — the  nowcast  analysis  domain  of  500  by  500  km 
is  located  at  the  center. 


Figure  3.  A  portion  of  the  upper-air  radiosonde  sounding  data  file,  new_sound. 


NOTE:  Columns  of  date  (from  left  to  right):  pressure  (mbar),  height 
above  sea  level  (m),  temperature  (°K),  dew-point  temperature  (°K), 
wind  speed  (m/s),  and  wind  direction  (deg). 
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2.2  Methodology 

2.2.1  Three-Dimensional  Analysis  of  MM5  and  Upper-Air  Data 

The  MM5  forecast  data,  as  shown  in  figure  1 ,  are  interpolated  to  the  three-dimensional 
nowcasting  grids  of  a  finer  spatial  horizontal  grid  resolution  than  that  of  the  MM5  and  to  the 
sigma-z  vertical  coordinate  used  by  the  current  IMETS  Battlescale  Forecast  Model  (BFM). 

This  vertical  coordinate  is  defined  as 


Z 


(1) 


Cartesian  vertical  coordinate, 
ground  elevation, 

material  surface  top  of  the  analysis  in  the  z*  coordinate,  and 
corresponding  height  in  the  z  coordinate  defined  by  H  =  H  +  zgmax 


maximum  value  of  the  terrain  elevation  in  the  nowcast  analysis  domain. 

Currently,  the  vertical  extent  of  the  nowcast  analysis  is  fixed  at  7000  m  above  the  highest  grid 
point  of  the  domain. 

The  following  procedures  are  taken  to  interpolate  the  raw  MM5  forecast  data  to  the  three- 
dimensional  grid  defined  as  the  nowcast  analysis  domain: 

•  Parameters,  including  temperature,  mixing  ratio,  and  horizontal  wind  components  (u  and  v 
on  the  pressure  level  pk),  are  vertically  interpolated  to  20  pre-determined  Cartesian-z 
vertical  levels  (zi,  i  =  1,20). 

•  Horizontal  interpolation  (to  analysis  i,  j  grid  space)  of  each  parameter  on  zi  levels  are 
performed  for  each  parameter. 

•  An  arbitrary  parameter  cp  on  zi  is  linearly  interpolated  to  a  z*  level  of  the  nowcasting  grid 
as 


where 


z  = 

the 

zg  = 

the 

H 

the 

H 

the 

where 

*Zj  g  max  — 

the 

<P 


where  ^  is  the  Cartesian  height  above  sea  level  of  z*,  calculated  from  eq  1  as 


(2) 
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(3) 


:Zg  +  Z 


■1//+Z 


gmax 


-zj 


H 


and  r  is  the  height  of  the  kth  z,  level. 

If  there  is  no  upper-air  radiosonde  sounding  data  available,  then  the  three-dimensional  data  of 

* 

(p^Z  )  determined  by  eq  2  is  used  as  three-dimensional  data  for  coupling  with  the  surface  field 
analyses. 

If  upper-air  radiosonde  sounding  data  are  available  at  the  time  of  analysis,  the  three-dimensional 
MM5  data  calculated  using  the  above  procedures  will  first  be  composed  with  the  upper-air 
radiosonde  sounding  data.  At  each  sounding  location,  radiosonde  data,  such  as  horizontal  wind 
vector  components,  temperature,  and  mixing  ratio,  are  vertically  interpolated  to  the  z,  levels  by 
using  a  linear  interpolation  method.  At  each  z;  level,  the  successive  correction  method  is  then 
performed  by  using  the  MM5  three-dimensional  data  as  background  data.  Finally,  linear  vertical 
interpolation  from  z;  to  zst  is  performed  for  each  parameter  and  grid  point.  At  this  point,  the 
composed  upper-air  univariate  analyses  are  ready  for  coupling  to  the  surface  field  univariate 
analyses  by  using  MARIAH  similarity  theory. 

2.2.2  Surface  Data  Analysis 

Using  the  MM5  forecast  as  background  data,  univariate  nowcasting  analyses  of  surface 
temperature,  dew-point  temperature,  and  horizontal  u  and  v  wind  vector  components  are 
conducted.  The  successive  correction  method  that  follows  is  used  in  this  development. 

If  cpa  (i,  j)  is  the  nowcast  value  at  the  grid  point  (i,  j),  and  cp  ;,(i,  j)  is  the  value  of  the  MM5 
background  at  the  same  grid  point,  then 


m 

<P.0’» = <PSUi) + 2  / - <pj 

k= 1 


(4) 


where 

m  is  the  observation  at  the  kth  location, 

/  o,k 

is  the  weight  for  each  observation, 

(O  is  the  value  of  the  background  at  the  observation  point  derived  by  a  bilinear  interpolation 

/  b,k 

method,  and 

m  is  the  total  number  of  observations. 
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The  weighting  factor  is  defined  as 


1 

Wk,ij  =  — 

rk,ij 


(5) 


where  yk  is  the  horizontal  distance  between  the  kth  observation  point  and  grid  point  (i,  j). 

The  method  is  repeated  in  an  iterative  fashion,  so  that  the  background  field  is  updated  by  the 
latest  analysis  after  each  iteration 


m 

<py  * » = -  <p.f  <»» 

k=\ 


(6) 


where 

(p  (n)  is  the  value  of  the  analysis  at  the  grid  point  after  the  nth  iteration,  and 
( 0  ( n )  is  its  value  interpolated  for  the  kth  location  after  the  nth  iteration. 

/  a,k 

In  this  development,  the  iteration  was  repeated  three  times,  for  practical  purposes,  to  obtain  the 
final  values.  A  future  modification  that  could  be  explored  would  be  to  add  a  vertical  “distance” 
component  to  the  yk  term  used  in  the  weighting  to  account  for  elevation  differences  between 

grid  point  and  observation  terrain  elevations. 

2.2.3  Similarity  Formulae,  MARIAH 

The  equations  and  approach  for  determining  the  similarity  scaling  constants  for  wind, 
temperature  and  specific  humidity,  referred  to  as  MARIAH  (Rachele  et  al.,  1995),  are  used  to 
complete  the  coupling  between  the  surface  and  upper-air  meteorological  analyses  in  order  to 
produce  the  final  three-dimensional  analysis  dataset.  For  practical  purposes,  it  is  assumed  that 
the  lowest  layers  above  ground  for  the  three-dimensional  analysis  data  are  fixed  to  the  standard 
screening  heights  used  for  collecting  most  surface  meteorological  data  (2  and  10  magi).  Details 
of  MARIAH  can  be  found  in  Rachele  et  al.  (1995);  here,  the  approach  is  only  briefly  described. 

Potential  temperature  0,  specific  humidity  q,  wind  speed  magnitude  V,  and  wind  direction  are 
calculated  from  the  MM5  lowest  sigma-P  level  data. 

The  scaling  height  z*  is  calculated  as 


z2  -  Zj  _  Az 
ln(z2  /  z, )  A  In  z 


(7) 
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and 


AV  =  V(z2)-V(z1) 
Ad  =  0(z2)-d(zl) 
Aq  =  q(z2)-q(zl) 


(8) 


Here,  z\  is  the  height  of  the  surface  observation  and  Z2  is  the  height  of  the  level  above  z\.  For  6 
and  q,  z\  is  2  m,  and  for  V,  z\  is  10  m. 

Depending  on  the  stability  of  the  atmosphere  near  the  surface,  the  profiles  of  temperature, 
specific  humidity,  and  wind  speed  are  determined  in  one  of  three  ways. 

For  the  Unstable  case  (A  6 1  Az  <  0),  the  Monin-Ohukhov  scaling  length  L  is  calculated  as 


L_  ^(Af)2 

\A6  +  0.61#A^]gAlnz 

(9) 

^=a-i5|ri/4 

(10) 

^=a-i5|ri/2 

(11) 

v  = 

(j)m  ln(z2  /  z, ) 

(12) 

e  -  kA0 

(j)H  ln(z2  /  z ,  ) 

(13) 

kAq 

q*  = - 

(j)H  ln(z2/zj) 

(14) 

Here,  k  is  Karman  constant  (=  0.4). 

Wind  speed,  potential  temperature,  and  specific  humidity  at  the  height  z2  are  given  as 


V*d) 

V(z2)  =  V(zl)  +  — ^Az 

KZ 

(15) 

d(z2)  =  d(z,)+  f  Az 

kz 

(16) 

q(z2)  =  q(zl)+  ^  Az 

(17) 
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For  the  stable  case  (At?/Az>0), 


i_^(z2)(l  +  0.61^(z2))(AF)2  5_ 

[A#  +  0.61#Ag]gAz  Zz 

(18) 

and 

fim  =  =  1  +  5(Z2  1  L) 

(19) 

F(z2  )  =  F(z, )  +  A  [ln(£l)  +  5 

/t  Z  j  T/ 

(20) 

9(z, )  =  ff(z, )  +  ^[ln(— )  +  5  — 

k  zx  L 

(21) 

q(z2)  =  q(zl)+  [ln(  2)  +  5  2  ‘] 

k  Zj  L 

(22) 

For  the  neutral  case  ( A#/  Az  =  0), 

F(z2)  =  F(z,)  +  Aln(£l) 

k  zx 

(23) 

9(z2)  =  0(z, )  +  -7-  ln(— ) 
k  zx 

(24) 

q(z2)  =  ^(zi)+~ln— 

k  Z, 

(25) 

3.  Example  Over  the  Southeastern  United  States 


To  demonstrate  the  ability  of  the  described  nowcasting  analysis  method,  it  was  applied  to  a 
domain  located  in  the  southeastern  United  States.  Figure  4  shows  the  terrain  elevation  contours 
across  the  500-by-500  km  nowcast  analysis  domain  (10-km  resolution),  centered  at  34.44°N  and 
84.55°W.  The  MM5  forecasts  were  obtained  directly  from  AFWA. 
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Figure  4.  Terrain  contours  for  the  southeastern  (SE)  U.S.  500-by-500  km 
nowcasting  area,  centered  at  34.44°N  and  84.55°W. 


The  symbol  represents  the  15-km  resolution  MM5  grid-point  locations,  the  letter  “S” 
represents  the  locations  of  the  surface  meteorological  observations,  and  the  letter  “U”  represents 
the  location  of  the  upper-air  radiosonde  observations.  A  radiosonde  site  (Nashville,  TN)  is 
located  in  the  northwestern  part  of  the  domain,  while  the  other  site  (Peach  Tree,  GA)  is  located 
in  the  south-central  part  of  the  figure. 

The  data  used  in  this  case  application  were  obtained  for  1200  UTC  on  April  12,  2004. 

Figure  5A  shows  the  surface  wind  vector  distribution  forecasted  by  the  15-km  resolution  MM5. 
In  figure  5  A,  the  wind  vectors  at  every  other  grid  point  are  plotted  for  visual  clarity;  the 
maximum  wind  vector  is  4.1  m/s.  Figure  5B  plots  the  observed  surface  wind  vectors  for  the 
corresponding  time. 
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Figure  5  A.  Surface  wind  vector  fields  for  the  10-km  resolution 
SE  U.S.  domain,  interpolated  from  MM5  15-km 
resolution  forecast  data,  for  1200  UTC,  April  12, 
2004. 


Figure  5B.  Observed  SE  U.S.  surface  wind  vectors  at 
1200  UTC,  April  12,2004. 
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The  nowcasting  analysis  result,  generated  by  applying  the  successive  correction  method  to  the 
wind  fields  shown  in  figures  5A  and  5B,  is  presented  in  figure  5C.  The  maximum  wind  vector  in 
the  figure  5C  is  4.0  m/sec. 


Figure  5C.  Nowcast  analysis  of  the  SE  U.S.  domain  surface 
wind  field,  using  the  successive  correction  method 
and  the  wind  fields  shown  in  figures  5A  and  5B. 

Figure  6A  shows  the  surface  temperature  field  forecasted  by  the  15-km  resolution  MM5,  valid  at 
1200  UTC  on  April  12,  2004.  The  nowcast  calculation,  achieved  by  combining  the  successive 
correction  method  and  MARIAH  formulae  using  the  MM5  forecast  temperature  field  and  the 
observed  temperature  data,  resulted  in  the  nowcast  analysis  temperature  field  shown  in 
figure  6B. 
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Figure  6 A.  Surface  temperature  field  for  the  10-km  resolution 
SE  U.S.  nowcast  domain,  interpolated  from  MM5 
15-km  resolution  forecast  data  for  1200  UTC, 
April  12,  2004. 


Figure  6B.  Nowcast  analysis  of  surface  temperature,  obtained 
from  the  successive  corrections  combination  of 
MM5  forecast  data  and  surface  observation  data  for 
1200  UTC,  April  12,  2004. 
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This  example  uses  upper-air  radiosonde  sounding  data  observed  at  both  the  Peach  Tree  and 
Nashville  sites. 

Figures  7  A  shows  temperature  and  dew-point  temperature,  and  figure  7B  shows  wind  speed  and 
direction,  obtained  for  the  Peach  Tree  radiosonde  site.  The  three  colored  curves  represent 
profiles  generated  from  the  15-km  MM5  forecast  (red),  the  actual  Peach  Tree  radiosonde 
observation  (purple),  and  the  successive  corrections  nowcast  analysis  (green). 


Figure  7A.  Vertical  distributions  of  temperature  and 
dew-point  temperature  for  the  Peach  Tree 
radiosonde  site,  generated  from  MM5, 
observed,  and  nowcast  data,  for  1200  UTC, 
April  12,  2004. 


Figure  7B.  Same  as  figure  7B,  except  for  wind  speed 
and  direction. 
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Similar  figures  were  obtained  for  the  Nashville  radiosonde  site,  but  are  not  shown  in  this  report. 

Figures  5-7  clearly  show  that  this  method  can  bring  three-dimensional  fields  of  meteorological 
parameters  forecasted  by  the  MM5  closer  to  the  actual  values  reported  by  observations. 


4.  Statistical  Evaluation  of  the  Nowcast  Analysis  Method  Over  the  Utah 
Western  U.S.  Domain,  Using  MesoWest  Surface  Observation  Data 


In  this  section,  the  nowcast  analysis  method  (for  surface  fields  only)  is  compared  statistically 
with  observed  data  obtained  from  the  MesoWest  surface  network,  for  a  domain  over  the  western 
U.S.,  centered  on  a  location  inside  Utah.  In  addition,  statistical  evaluations  are  made  for  both  the 
pure  MM5  forecast  and  for  a  separate  method  based  on  the  earlier  BFM  dynamic  assimilation  (or 
nudging)  methodology. 

The  three  different  set  of  results  being  compared  in  this  section  are  as  follows: 

1.  Surface  meteorological  fields  generated  by  the  nowcast  analysis  method,  using  MM5, 
domain  2  (15-km  resolution),  12-h  forecasts  as  background  data 

2.  12-h  forecasted  surface  meteorological  fields  from  the  MM5,  domain  3  (5-km  resolution) 

3.  Surface  meteorological  fields  produced  by  the  BFM  nudging  method,  using  a  3-h  dynamic 
assimilation  cycle  (Henmi  and  Dumais,  1998) 

In  the  above  comparisons,  a  5-km  resolution  grid  was  used  to  produce  the  final  datasets. 

4.1  Modeling  Domains 

The  three  different  methods  were  applied  to  a  5-km  horizontal  resolution  domain  centered  on 
Utah  (fig.  8). 
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Figure  8.  Geographical  extent  of  the  three  MM5  computational  nests. 

The  MM5  model  was  executed  in  a  triple-nest  configuration.  The  model  domains  1  and  2  have  a 
horizontal  mesh  size  of  55-by-55  grid  points,  while  the  innermost  model  domain  3  has  a 
horizontal  mesh  size  of  85  by  85.  The  horizontal  resolutions  of  each  domain  are  45,  15,  and  5 
km,  respectively.  For  this  study,  the  forecast  fields  for  domain  2  were  used  as  the  background 
data  for  the  nowcast  analysis  algorithm  (and  for  the  BFM  nudging  method),  since  1 5  km  is  the 
smallest  resolution  of  MM5  data  that  is  typically  available  in  IMETS. 

Figure  9  shows  the  terrain  contours  across  the  grid  domain  used  for  the  statistical  comparisons, 
which  has  5  l-by-5 1  grid  points  and  a  5-km  horizontal  resolution.  The  letter  “S”  represents  the 
location  of  MesoWest  surface  observation  sites,  whose  data  were  used  for  producing  the  actual 
nowcast  (and  BFM  nudging-based)  analyses,  and  the  symbol  represents  those  observation 
sites,  whose  data  were  used  only  to  compare  to  the  nowcast  (or  forecast/BFM  nudging)  analyses 
to  produce  the  statistics  (i.e.,  independent  verification  observations). 
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Figure  9.  Terrain  contours  for  the  5-km  resolution  Utah  study 
domain  and  the  locations  of  MesoWest  surface  data. 


4.2  Meteorological  Data 

As  shown  in  figure  9,  numerous  MesoWest  surface  observational  data  sites  are  available.  The 
data  are  publicly  available  from  the  University  of  Utah  MesoWest  cooperative  anonymous  file 
transfer  protocol  (FTP)  site.  MesoWest  data,  which  is  archived  daily  at  ARL  into  an  ASCII  file 
called  total .  dat,  includes  observations  taken  over  the  past  24-h  period,  in  15-min  intervals, 
for  the  entire  western  United  States.  A  detailed  description  of  the  MesoWest  data  can  be  found 
in  Horel  et  al.  (2002). 

The  National  Center  for  Environmental  Prediction  (NCEP)  Aviation  Model  (AVN)  provides 
global  gridded  forecast  information  at  a  1°  resolution.  These  data  were  used  for  initialization  and 
time-dependent  lateral  boundary  conditions  for  the  locally  run  MM5  model.  The  AVN  forecast 
data  initialized  at  1800  UTC  were  used  for  this  study  and  were  obtained  through  the  anonymous 
FTP  site  at  NCEP.  The  use  of  AVN  (now  called  GFS)  is  consistent  with  the  external  model  used 
by  AFWA  to  provide  lateral  boundary  conditions  for  its  operational  MM5  runs. 

Datasets  were  collected  for  two  different  study  periods,  during  which  successive  correction 
nowcast  analysis,  BFM  nudging  assimilation,  and  MM5  forecast  data  were  generated  and 
statistically  analyzed  for  the  following  periods: 

•  August  1  through  21,  2003,  for  17  days 

•  December  31,  2003  to  January  28,  2004,  for  17  days 
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MM5  15-km  resolution  forecast  data,  as  shown  in  figure  1  and  used  for  analysis  background 
data,  are  stored  at  output  frequency  of  every  three  hours  (0000,  0300,  0600,  0900,  1200,  1500, 
1800,  and  2100  UTC).  Similarly,  sf  cmet .  d  surface  data  files,  as  shown  in  figure  2,  are  made 
every  three  hours,  corresponding  to  the  MM5  output  hours  listed  above. 

4.3  Statistical  Parameters 

The  following  statistical  parameters  were  calculated  for  evaluating  the  quality  of  the  successive 
corrections  nowcast  analysis,  BFM  nudging  assimilation,  and  MM5  5 -km  resolution,  12-h 
forecast  datasets: 

•  Mean  difference  (MD) 

•  Mean  absolute  difference  (AD) 

•  Root  mean  square  error  (RMSE) 

•  Root  mean  square  vector  error  (RMSVE) 

•  Correlation  coefficient  (CC) 

•  Mean  wind  direction  difference  (MWDDF) 

Details  of  these  parameters  are  described  in  Henmi  (2003). 

For  this  report,  these  statistical  parameters  were  calculated  separately  for  the  two  different  study 
periods,  inclusive  of  all  17  days  and  all  analysis  periods  (8)  per  day. 

4.4  Evaluation  Results 

Statistical  parameters  were  calculated  between  the  surface  nowcast  analysis  (and  MM5 
forecast/BFM  nudged)  data  and  surface  observation  data  (after  interpolation  to  the  5-km 
resolution  study  domain)  and  were  obtained  at  the  sites  marked  using  the  symbol  as  shown 
in  figure  9.  Tables  1-3  show  the  values  of  statistical  parameters  calculated  for  the  period  in 
August  2003. 

Statistical  parameters  calculated  for  the  period  in  January  2004  are  shown  in  tables  4-6. 
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Table  1.  Statistical  parameters  between  the  successive  corrections  nowcast  surface  analysis  (using  MM5 
15 -km  resolution,  12-h  forecast  background  fields)  and  surface  observation  data  for  August  2003 
data. 


MD 

AD 

RMSE 

cc 

Temperature  (°C) 

-1.0 

3.0 

3.9 

.82 

Dew-point  temperature  (°C) 

-0.9 

2.3 

3.2 

.75 

Wind  speed  (m/sec) 

-0.9 

1.9 

2.5 

.36 

Wind  vector  x-component  (m/sec) 

1.1 

2.1 

3.0 

.13 

Wind  vector  y-component  (m/sec) 

-0.5 

2.1 

2.7 

.57 

NOTE:  RMSVE  =  4.0  m/sec,  MWDDF  =  49° 

Table  2.  Statistical  parameters  between  the  12-h  forecast  surface  data  from  the  5-km  resolution  MM5, 
domain  3,  and  surface  observed  data  for  August  2003  data. 


MD 

AD 

RMSE 

cc 

Temperature  (°C) 

-1.7 

3.1 

3.9 

.81 

Dew-point  temperature  (°C) 

4.2 

5.2 

7.1 

.27 

Wind  speed  (m/sec) 

0.8 

2.1 

2.7 

.36 

Wind  vector  x-component  (m/sec) 

0.7 

2.3 

3.3 

.03 

Wind  vector  y-component  (m/sec) 

1.4 

2.6 

3.3 

.52 

NOTE:  RMSVE  =  4.6  m/sec,  MWDDF  =  53° 


Table  3.  Statistical  parameters  between  the  BFM  nudging  assimilated  surface  fields  (using  MM5  15-km 

resolution,  12-h  forecast  background  fields)  and  surface  observed  data  for  August  2003  data.  The 
BFM  assimilation  period  was  3  h. 


MD 

AD 

RMSE 

cc 

Temperature  (°C) 

-2.3 

3.3 

3.9 

.85 

Dew-point  temperature  (°C) 

3.8 

4.3 

5.4 

.60 

Wind  speed  (m/sec) 

-1.1 

1.9 

2.5 

.44 

Wind  vector  x-component  (m/sec) 

0.0 

1.8 

2.5 

.20 

Wind  vector  y-component  (m/sec) 

0.1 

2.1 

2.8 

.48 

NOTE:  RMSVE  =  3.7  m/sec,  MWDDF  =  54° 
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Table  4.  Same  as  table  1,  except  for  January  2004  data. 


MD 

AD 

RMSE 

CC 

Temperature  (°C) 

-2.0 

3.7 

6.9 

.57 

Dew  point  temperature  (°C) 

-0.4 

2.5 

6.2 

.60 

Wind  speed  (m/sec) 

-0.4 

1.9 

2.9 

.35 

Wind  vector  x-component  (m/sec) 

0.5 

1.7 

2.6 

.28 

Wind  vector  y-component  (m/sec) 

-0.8 

2.0 

2.9 

0.40 

NOTE:  RMSVE  =  3.9  m/sec,  MWDDF  =  55° 


Table  5.  Same  as  table  2  except  for  January  2004  data. 


MD 

AD 

RMSE 

CC 

Temperature  (°C) 

0.8 

3.6 

6.9 

0.50 

Dew  point  temperature  (°C) 

0.5 

4.9 

8.2 

0.39 

Wind  speed  (m/sec) 

1.6 

2.6 

3.4 

0.35 

Wind  vector  x-component  (m/sec) 

1.1 

2.0 

3.2 

0.18 

Wind  vector  y-component  (m/sec) 

0.4 

2.4 

3.3 

0.50 

NOTE:  RMSVE  =  4.5  m/sec  MWDDF  =  42° 

Table  6.  Same  as  table  3  except  for  January  2004  data. 


MD 

AD 

RMSE 

CC 

Temperature  (C) 

0.8 

3.1 

6.4 

0.55 

Dew  point  temperature  (  C) 

-0.1 

2.9 

6.4 

0.55 

Wind  speed  (m/sec) 

-0.2 

1.8 

2.6 

0.46 

Wind  vector  x-component  (m/sec) 

-0.1 

1.6 

2.5 

0.30 

Wind  vector  y-component  (m/sec) 

0.4 

1.9 

2.7 

0.48 

NOTE:  RMSVE  =  3.6  m/sec,  MWDDF  =  43° 
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The  statistical  results  shown  in  tables  1-6  were  obtained  by  combining  the  data  over  all  eight 
daily  analysis  periods  (0000,  0300,  0600,  0900,  1200,  1500,  1800,  and  2100  UTC. 

Figures  10  and  1 1  are  scatter  diagrams  for  temperature — the  first  compares  the  successive 
corrections  nowcast  analysis  results  with  observation  data,  and  the  second  compares  the  12-h 
forecast  data  from  the  MM5,  domain  3  (5-km  resolution),  with  observation  data.  All  diagrams 
are  for  surface  data  only  and  for  results  on  the  Utah  5-km  resolution  grid  mesh.  Similarly, 
figures  12  and  13  are  the  scatter  diagrams  for  dew-point  temperature,  and  figures  14  and  15  are 
the  scatter  diagrams  for  wind  speed.  These  figures  are  for  the  period  in  August  2003. 


Figure  10.  Scatter  diagram  between  successive  corrections  nowcast 
analysis  surface  temperature  data  with  and  surface 
temperature  observations,  for  the  period  in  August  2003 
over  Utah  domain. 
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Figure  1 1 .  Scatter  diagram  between  pure  12-h  MM5,  5-km  resolution 
forecast  surface  temperature  data  and  surface  temperature 
observations,  for  the  period  in  August  2003  over  Utah 
domain. 
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Comparisons  of  tables  1-6  and  figures  10-15  show  that  meteorological  fields  results  can  be 
improved  by  applying  the  successive  corrections  nowcast  analysis  method,  beyond  the  results 
that  an  earlier-run,  high-resolution  numerical  weather  prediction  (NWP)  model  can  produce. 
However,  over  extremely  complex  terrains,  such  as  the  Utah  study  domain,  the  current  method 
indicates  only  a  limited  potential  for  vastly  improving  surface  meteorological  fields  (especially 
winds).  In  a  similar  study  over  Oklahoma,  where  the  terrain  features  were  not  as  complex 
(Sauter  et  al.,  2001),  results  using  this  method  produced  far  improved  surface  meteorological 
fields  over  the  coincident  NWP  model  forecast  calculations. 

The  following  time  series  of  the  statistical  parameters  (figs.  16-23)  show  MD,  CC,  AD, 
RMSVE,  and  MWDDF,  calculated  for  every  three  hours  of  the  diurnal  cycle,  using  all  available 
data  for  the  5-km  resolution  Utah  domain.  Figures  16-19  are  for  the  period  in  August  2003,  and 
figures  20-23  are  for  the  period  in  January  2004.  In  these  figures,  the  statistical  parameters 
calculated  for  the  three  different  methods  (successive  corrections  nowcast  analysis;  5-km 
resolution,  12-h  MM5  forecasts;  and  BFM  nudging  assimilation  over  a  3-h  period)  are 
distinguished  using  three  different  colored  lines. 


Figure  16.  Time  series  of  statistical  parameters  for  all  methods  (MD,  CC,  and  AD),  for  surface 
temperature,  for  the  period  in  August  2003  in  the  Utah  domain. 
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Figure  17.  Same  as  figure  16,  except  for  surface  dew-point  temperature. 


Figure  18.  Same  as  figure  16,  except  for  surface  wind  speed. 
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Figure  19.  Same  as  figure  18,  except  for  RMSVE  and  MWDDF. 


Figure  20.  Time  series  of  the  statistical  parameters  for  all  methods  (MD,  CC,  and  AD),  for 
surface  temperature,  for  the  period  in  January  2004  in  the  Utah  domain. 
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Figure  22.  Same  as  figure  20,  except  for  surface  wind  speed. 
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Figure  23.  Same  as  figure  22,  except  for  RMSVE  and  MWDDF. 

As  seen  from  these  figures,  both  the  successive  correction  nowcast  analysis  and  the  BFM 
nudging  assimilation  methods  (which  incorporate  recent  local  observations  into  their  solutions) 
clearly  produced  better  statistical  results  than  the  pure  12-h  forecasts  obtained  from  the  MM5, 
domain  3  (5-km  resolution).  The  differences  between  the  nowcast  analysis  and  BFM  nudging 
methods  appear  to  be  small  statistically.  It  is  worth  noting  that,  during  both  periods,  the 
statistical  parameters  obtained  by  all  three  methods  were  not  significantly  high.  As  seen  in  the 
statistical  parameters  of  forecast  results  by  MM5,  domain  3,  the  MM5  did  not  perform  well  over 
the  complex  domain  of  Utah  for  either  the  warm  (August  2003)  or  the  cold  (January  2004) 
period. 


5.  Concluding  Remarks 


This  report  focuses  on  a  nowcast  analysis  method,  based  on  successive  corrections  and 
developed  to  use  MM5  forecast  fields  as  background  data  to  combine  surface  and  upper-air 
radiosonde  observation  data  (Henmi,  2003),  which  has  been  modified  for  easier  integration  into 
the  U.S.  Army  IMETS.  This  method  utilizes  MM5  GRIB  forecast  data,  surface  observations, 
and  upper-air  radiosonde  observations  similar  to  those  currently  obtained  from  the  AFWA  in 
IMETS.  This  method  starts  by  composing  and  univariately  analyzing  two  different  sets  of  data: 
first,  the  surface  observations  and  the  corresponding  MM5  surface  forecast  data,  then  the  upper- 
air  radiosonde  sounding  data  and  the  MM5  forecast  data.  After  these  two  steps,  the  similarity 
theory  formulae,  MARIAH,  are  used  to  couple  the  surface  data  univariate  analysis  fields  with  the 
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corresponding  upper-air  analysis  fields,  by  calculating  the  meteorological  parameters  at  the 
levels  near  the  ground  using  the  MARI  AH  formulae. 

The  usefulness  of  the  method  was  demonstrated  by  applying  it  to  a  domain  centered  in  the 
southeastern  United  States,  for  which  MM5  GRIB  forecast  data,  surface  observations,  and  upper- 
air  radiosonde  observations  were  available  from  AFWA. 

Evaluation  and  comparison  studies  of  three  different  methods  were  conducted  over  a  250-by- 
250  km  domain  in  Utah,  and  the  results  of  each  method  were  statistically  comparing  with  the 
Utah  MesoWest  surface  data.  The  three  different  methods  compared  were  as  follows: 

•  Successive  corrections  nowcast  analysis  data,  which  used  the  MM5,  domain  2  (15-km 
resolution),  12-h  forecasts  as  background  data 

•  Pure  12-h  forecast  data  produced  by  the  MM5  domain  3  (5-km  resolution) 

•  BFM  nudging  assimilation  method,  using  a  3-h  assimilation  period 

Two  observation  datasets,  for  August  2003  and  January  2004,  were  statistically  compared  with 
the  above  three  methods. 

The  two  non- forecast  methods  produced  superior  statistical  results  as  compared  to  the  pure 
forecast  method;  however,  there  were  no  significant  statistical  differences  between  the  two  non- 
forecast  method  datasets  themselves.  Statistical  results  showed  that  nowcasting  of  the 
temperature  field  can  be  accomplished  more  accurately  than  nowcasting  of  the  dew-point 
temperature  and  wind  fields.  Over  a  complex  terrain,  such  as  in  the  model  domain  used  for  this 
study,  forecasts  and  nowcasts  of  wind  fields  proved  to  be  difficult,  probably  due  to 
inhomogeneous  distribution  of  observation  sites  and  the  limitations  of  the  forecast  model.  It  is 
essential  to  have  better  forecast  data  as  background  data  in  order  to  use  the  nowcast  analysis 
method  more  effectively. 
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Acronyms 


AD 

Mean  Absolute  Difference 

AVN 

Aviation  Model 

BFM 

Battlescale  Forecast  Model 

CC 

correlation  coefficient 

FTP 

file  transfer  protocol 

MD 

mean  difference 

MM5 

Penn-State/NCAR  Mesoscale  Model  version  5 

MWDDF 

mean  wind  direction  difference 

NCAR 

National  Center  for  Atmospheric  Research 

RMSE 

root  mean  square  error 

RMSVE 

root  mean  square  vector  error 

SE 

southeastern 
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